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SECOND QUARTERLY PROGRESS REPORT UNDER CONTRACT NAS5-27334 


STUDY OF LANDSAT-D THEMATIC MAPPER PERFORMANCE 
AS APPLIED TO HYDROCARBON EXPLORATION 


Period Covered: January 7 , 1983 to April 7, 1983 



During this reporting period, EarthSat was able to perform a first sub- 
stantial analysis of Thematic Mapper data, and to make a first evaluation 
of its potential use for hydrocarbon exploration. The results of this 
work, together with analytical work concerning band inage selection for 
Thematic Mapper data, are contained in three reports: (1) "Evaluation of 

Haematic Mapper Performance as Applied to Hydrocarbon Exploration, " 

(2) "Inpli cations of Information from Lands at- 4 for Private Industry," 
and (3) "Selecting Band Combinations with Thematic Mapper Data" (attached) . 
The following summary canxoents are in large part taken from the above 
reports. Presentation of these results was made at the Landsat Early 
Results Symposium (February 22-24, 1983) and the Goddard Memorial Symposium 
(March 24-25, 1983) . 


II. A. Analysis of a July 25, 1982 scene of southern Ontario shews that 
TM data permits improved delineation of known oil and gas fields. 
Visible lineaments appear to be the surface expression of major 
through-going structures, controlled by subsurface fracturing, not 
readily visible on earlier MSS data. 

B. The use of hue, saturation, and value (HSV) image processing 
techniques on the Death Valley scene of 17 November 1982 permits 
direct comparison of 134 processed imagery with existing 1:250,000 
scale geological maps of the area. In particular, small outcrops 
of Tertiary volcanic material overlying Paleozoic sections shew a 
detail not provided by the geologic maps. ‘This type of comparative 
mapping is not possible using natural color or standard false color 
presentations of TM data. In fact, the processing employed all 

7 TM bands, with a TM 5/TM 2 ratio controlling hue, a HI 5/TM 7 
ratio controlling saturation, and the first principal component 
combination controlling the value. 

TM 2 was doc sen for its sensitivity to iron oxides, TM 7 (2.2 
microns) for its sensitivity to hydroxyl bands, and TM 5 for its 
high variance. 

C. Analysis has just begun of TM data over Lawton, Oklahoma. We 
believe that we see evidence that the reducing diemical environment 
associated with hydrocarbon microseepage has changed ferric iron 
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to soluble ferrous iron, allowing it to be leached from the 
Rush Springs Sandstone, with a consequent bleaching of the sur- 
face. To date, this work has been done using only interactive 
computer displays, and the next stage of analysis will be the 
production of enhanced hard copy outputs. 

D. Hie Owl Creek, Wyoming scene has been received, but the light 
dusting of snow that covers much of the area will, we believe, 
hinder spectral analysis of this imagery. However, the altera- 
tion on this scene associated with the Copper Mountain uranium 
deposit is visible. The overall configuration of the 0/1 Creeks 
suggest that they lie within a major left-lateral wrench zone 
broken by northwest trending thrust faults paralleling the Wind 
River thrust to the southwest. We regard the amount of structural 
information available in this scene as spectacularly high. 

E. The band selection algorithm referenced above lias been applied 
to each scene and each subscene analyzed. The results show a 
surprising consistency, with the 1,4,5 combination selected as 
optimal in almost all cases. This has been true in scenes as 
diverse as Washington, D.C. , Death Valley, and Lawton, Oklahoma. 
Hie sane computation shows that both the natural color (1,2,3) 
and the standard false color (2,3,4) combinations lie far down in 
the rankings. The production of 1,4,5 hard copy images confirms 
that this combination usually provides more and more easily 
interpretable information than either the natural or standard 
fait -* color images, at least so far as geological analysis is con- 
cerned. 


III. General Comments and Problem Areas 


Discussions at the February Early Results Symposium make it clear that many 
investigators are exploring the radiometric and geometric properties of TM 
imagery and evaluating the performance of the sensor system, but no other 
group is evaluating the instrument's performance in the context of hydro- 
carbon exploration. We therefore think it logical tliat we should enphasize 
the latter aspect, while working closely with other investigators to minimize 
duplication in engineering evaluation of the system. 

In this context, two main factors are presently limiting our actixities . 
First, we have found it difficult to obtain cloud-free, snow-free imagery 
over the test site areas of interest for hydrocarbon exploration. Second, 
our analysis to date lias been limited to P tapes. Since a number of our 
evaluation activities depend upon radiometric modification of line data, 
access to A tapes is important. 

Continued analysis of test site data for structural features, evidence of 
hydrocarbon microseepage, and botanical anomalies, together with processing 
of A tape data, will constitute major activities of the next quarterly 
period. 


EVALUATION OF THEMATIC MAPPER PERFORMANCE 
AS APPLIED TO HYDROCARBON EXPLORATION 


Dr. John R. Everett, Dr. Charles Sheffield, Dr. Jon Dykstra 
Earth Satellite Corp., 7222 47th Street, Chevy Chase, MD 20815 


Results from a limited sample of Landsat-4 TM scenes demonstrate that the 
greater number of spectral bands and the increased spectral and spatial re- 
solution over previous satellite data greatly expand, improve and refine the 
geological inferences possible from space. The number and narrower width of 
available spectral bands, coupled with digital processing techniques such as 
band ratioing, principal component analysis and hue-saturation-intensity pre- 
sentation, allow the differentiation and mapping of a large number of geologic 
units and the identification of altered rock associated with intrusions in the 
Death Valley, California area. These capabilities will contribute greatly to 
geologic exploration in arid areas. 

The increased spatial resolution also assists in distinguishing glacial from 
tectonic features (difficult on previous imagery) in the Detroit area. Several 
of the Ordovician oil and gas fields of Ontario lie on fracture zones clearly 
visible on Landsat-4 imagery. All of these fields rely on fracture-related 
porosity and permeability for production. The current-borne sediment plumes 
in Lake St. Clair and Erie demonstrate that Landsat-4 will contribute to our 
understanding of lacustrine, riverine and marine depositional and erosional 
patterns and processes. 

Digitally enhanced Landsat-4 TM imagery reveals geological features of the 
Anadarko Basin, Oklahoma, in greater numbers and detail than previously 
possible. Comparison with existing work suggests that some of the subtle 
spectral differences detected are related either to changes in chemical compo- 
sition or resultant changes in vegetation related to micro-seepage of hydro- 
carbons adjacent to oil fields in the area. 

There seem to be several categories of features (throughgoing structures, 
large area subtle tonal contrasts, etc.) that are still best seen on older 


Landsat imagery. Thus, Landsat-4 does not make earlier systems obsolete; 
rather, it provides a powerful new tool to exploration geologists. 

There is a fortuitous synchronous development of satellite remote sensing, 
and the concepts of plate tectonics and vertical migration of hydrocarbons - 
two paradigms which are reshaping geological exploration thinking. The 
steadily improving quality of satellite data is providing geologists with a 
valuable tool with which to test, revise and exploit these new concepts. 


EVALUATION OF THEMATIC MAPPER PERFORMANCE 
AS APPLIED TO HYDROCARBON EXPLORATION 


Dr. John R. Everett, Dr. Charles Sheffield, Dr. Jon Dykstra 
Earth Satellite Corp., 7222 47th Street, Chevy Chase, MD 20815 


EVALUATION OF THEMATIC MAPPER PERFORMANCE 
AS APPLIED TO HYDROCARBON EXPLORATION 


TM AND GEOLOGIC EXPLORATION 

Since the early 1 960 ' s , the science of geology has been undergoing a major 
revolution. The new paradigm of plate tectonics and seafloor spreading is 
replacing the older paradigm of a rigid, stable earth. Inherent in the 
acceptance of plate tectonic theory is a growing appreciation of the role of 
plate motion in determining the location of mineral deposits and hydrocarbon 
accumulations. It is fortunate that developments in spaceborne remote 
sensing have paralleled these developments in geologic thinking. As a con- 
sequence, we have remote sensing tools that view the earth with appropriate 
scale and scope to enable us to appreciate and map the regional structures 
that reflect the motions of continent-sized segments of the earth's crust. 

We received our first glimpses of the earth from space with photos from the 
Apollo and Gemini flights. The first three Landsat satellites gave us near 
ubiquitous high resolution (80 metre) coverage in four spectral bands. These 
data have had and continue to have enormous impact on all facets of the per- 
ception and management of renewable and non-renewable natural resources and 
the environment. 

Before speculating on the impact of the new data types from Landsat-4, it is 
useful to take a look at the role data from the first three Landsat satellites 
have in geologic exploration and their current level of acceptance. In a 
general sense, Landsat data has made its major contribution to hydrocarbon 
exploration in the spatial domain. In mineral exploration, Landsat has 
revealed some spectral information, but again the major contribution is 
spatial. The synoptic view of over 34,000 square kilometres of the earth's 
surface on a single Landsat image permits the detection and mapping of major 
regional structures associated with the geologic development of entire geo- 
logic provinces. It is also possible, through special digital enhancements, 
to map some of the more subtle surface expressions of fracturing, folding and 
alteration associated with buried hydrocarbon accumulations and the emplacement 


of mineral deposits. The data make it possible to interrelate widely sepa- 
rated geologic features and detect subtle changes that occur over tens of 
miles and, hence, have gone unnoticed on conventional types of data. Perhaps 
most important of all, the new perspective that the view from space provides 
stimulates us, even forces us, to think of geology in new ways and perceive 
new possibilities. Truly it isn't a panacea but it is an extremely powerful 
tool. It has not by itself "found" an oil field or mine, but it has made 
significant contributions to the exploration thinking that led to the discovery 
of millions of barrels of oil and millions of tons of ore. 

At this point, it is appropriate to consider how industry regards this new 
tool. Bob Porter (President, Earth Satellite Corporation) once characterized 
the oil and gas and mineral exploration industries' acceptance of Landsat 
data or any technical innovation as progressing through three stages: 

1. An initial "Gee Whiz" stage; 

2. An "Interesting but not yet commercial" stage; and finally 

3. Full acceptance and integration of the innovation as a bona fide 
exploration tool. 

The progression through these stages appears to be a direct function of the 
geologists' and geophysicists' experience and familiarity with the tool; in 
this instance, satellite data. Consequently, not all exploration groups 
progress through each stage in unison. 

In the earlier "Gee Whiz" stage, the geologist may be impressed by the 
presence of a familiar structure or landform that is visible on the imagery, 
or by a subtle feature visible on an exotic digitally enhanced color image. 

At this stage, some are tempted to see satellite data as the answer to all 
exploration woes, and the potential to "oversell" its capabilities is a real 
danger. 

This stage is usually quickly replaced by the "Interesting but not yet 
commercial" stage where imagery produced from the Landsat data is viewed 
mainly as a cheap, low resolution, aerial photo substitute. To so view 
Landsat imagery is to miss both the unique geologic perspective of the 
imagery and the potential contribution available through computer processing 
and data base integration of the digital data. 

The final stage is characterized by the integration of the geologic information 
gleaned from the Landsat imagery with a complete oil and gas or minerals 
exploration model. At this point, the exploration geologist fully appreciates 
the potentials and limitations provided by the satellite system and routinely 
applies the satellite data to the design and execution of exploration pro- 
grams. This whole process takes time. It took the gravity meter and the 
seismograph about 40 years to make the transit; the digital computer has 
progressed more rapidly. 

At present, I think a healthy percentage of the exploration community is 
working out of the second stage of acceptance. There is a great promise that 
the Thematic Mapper data will provide the results necessary to convince the 
exploration geologist and, more importantly, the exploration financial managers. 


to move fully into the final stage. However, in order for this to occur, 
there must exist a continuity of data and, equally important, an availability 
of data at a reasonable (that is, user justifiable) price. 

The two major advantages of Thematic Mapper data over that of the MSS system 
are its increased spatial resolution and its greater number of narrow, 
strategically placed, spectral bands (Table 1). The 30 metre pixel size will 
permit finer definition of ground features and thereby improve the reliability 
of photo-geologic interpretation of geologic structure. Of equal importance 
is the increased homogeneity of the types of surface material within a given 
pixel. The less mixed the pixel, the greater the potential of extracting 
useful spectral information. The increased spectral resolution is allowing 
geologists to map altered zones associated with mineralization based not only 
on iron oxides, but on the basis of recognizing rocks and soils rich in 
hydroxyl groups, such as many of the clays formed as a product of the minerali- 
zation process. 

The increased spectral sensitivity also promises the ability to detect some 
types of vegetation changes that are associated with anomalous mineralization. 
This will be particularly helpful where soil and plants obscure the bedrock. 

This capability is not definitely proven, but it is theoretically possible 
and highly anticipated. 

In addition to plate tectonics, there is a second revolution going on in the 
geologic thinking of petroleum exploration. The old paradigm of tightly 
sealed hydrocarbon traps which retain for long periods of time petroleum that 
was generated and migrated in the distant past is giving way to a newly 
evolving paradigm which envisions a much more dynamic scenario in which most, 
if not all, traps leak, and the generation and migration of hydrocarbons is a 
continuing process. This Implies that there is very little, if any, really 
old oil or gas, rather, only new hydrocarbons generated from old rocks or 
retained in old traps. The hydrocarbon leaked from these imperfect traps 
moves vertical! ly through the overlying rocks to the surface and, in the 
course of its movement, produces a host of chemical changes. The near surface 
environment manifests this leakage in a variety of geochemical, biological, 
geobotanical , or geomorphological anomalies and by the simple presence of 
hydrocarbon itself. 

This new paradigm also has Important significance to the mineral explorationists. 
The chemical environment created by leaking hydrocarbon has caused the emplace- 
ment of a vast amount of lead, zinc, uranium, and silver and has potentially 
played a role in localizing some deposits of gold, copper, and barite. 

Iri the remainder of the talk, I would like to give you an appreciation of a 
few of the ways that the increased spectral and spatial resolution of the 
Thematic Mapper will affect geologic exploration. First, the increased 
spatial resolution. It is clear to us that, with careful computer and photo- 
graphic processing, the quality of the TM digital data enables sharp photo- 
graphic enlargements to a scale of 1:50,000 and, in some cases, larger. With 
clearly interpretable Imagery at these scales, exploration geologists are 
able to significantly refine their structural interpretations compared to 
those made from 80 metre resolution MSS imagery. 


0 What was detectable only as a lineament on MSS imagery might be 
able to be confidently mapped as a fault. 

0 More importantly, some of the smaller features that indicate 

direction of movement along faults will be Identifiable in the TM 
imagery, where they are lost in the lower resolution MSS imagery. 

Southern Ontario (Example) 


There are relatively few TM scenes available of areas with strong oil and gas 
interests, however, the July 25 image of Detroit does include some of the oil 
and gas fields of southern Ontario, Canada. 

Thematic Mapper False Color Composite of Southern Ontario 

On the false color image of the Ontario area, we have delineated some of the 
more prominent 1 inears along with the location of the Malden, Colchester and 
Leamington oil and gas fields. At the Malden and Colchester fields, the 
hydrocarbon accumulations are in fractured, dolomitized, Ordovician limestones. 
The fractures trend WNW. It's a safe bet that the lineament marked on the 
imagery is the surface expression of a major through-going structure which is 
controlling the subsurface fracturing. 

The Leamington field is a little younger and is located in ancient reef 
deposits. ^Reefs are known to prefer the high edge of structural blocks. The 
intersecting lineaments mapped on the imagery may well define two intersecting 
normal fault zones responsible for the uplift of a block edge and the locali- 
zation of the Silurian reef. 

Cement, Oklahoma (Example) 

The most exciting potential contribution of the TM data is the availability 
of seven carefully placed spectral bands. For oil and gas exploration, these 
spectral bands will be extremely useful for the detection of possible surface 
rock alteration and geobotanical anomalies associated with microseepage of 
hydrocarbons from buried oil and gas accumulations. As the leaking hydro- 
carbons find their way to the surface, they alter the chemistry of the rocks 
through which they pass. At the surface, several things can occur: the 

altered surface chemistry may change the spectral and erosional characteristics 
of the surface rocks and soils and/or it may create variations in the density, 
type or vigor of any vegetation growing in the altered soils. 

As a portion of our TM investigations, we are assessing TM's ability to 
identify and deli reate areas of surface alteration due to microseepage of 
hydrocarbons at the Cement oil and gas field. The Cement field is overlain 
on the surface by the Rush Springs Sandstone. The Rush Springs Sandstone is 
a characteristic red color due to an abundance of ferric iron-oxides. The 
reducing chemical environment associated with the hydrocarbon microseepage 
appears to have altered the insoluble ferric iron to soluble ferrous iron, 
allowing it to be leached out of the sandstone. The result is a bleaching of 
the sandstone overlying the oil and gas field. We have only recently received 
our TM coverage of the Cement area. Our spectral analyses are therefore in 
their earliest stages. These next few images are pictures of our interactive 


Image processing system and demonstrate our first cut approach at delineating 
the bleached area associated with the Cement field. 

Cement NCC/HSV 

The iron oxide rich areas are delineated of the imagery as areas of orange to 
red color. The quarter frame image shows the location of the red coloration 
following closely the outcrop of the Rush Springs Sandstone. However, the 
area overlying the Cement oil and gas field is one of those areas where the 
Rush Springs appears to have lost its strong iron oxide signature. 

DEATH VALLEY, CALIFORNIA (EXAMPLE) 

The first place to assess potential of TM data to map different rock types is 
in an area of low vegetation density and diverse rock types. Within the 
present range of the TM system there is clearly no more vegetation free area 
than Death Valley, California. 

(Death Valley NCC/Death Valley FCC) - (Eigen/HSV) 

The following imagery are of an approximately 1/3 TM scene area of the Death 
Valley, California overpass on 17 November, 1982. The scene includes a 
natural color, false color, eigen and HSV image. 

The Hue, Saturation and Value (HSV) image is one of the more exciting images 
for geologic applications. Through the use of two ratios as hue and saturation, 
and the first eigenband as the value, the resulting HSV image possesses the 
spectral information of a ratio image and the spatial integrity of the first 
eigenband. 

The hue of the image is controlled by the ratio of TM5 (1.6 microns) over TM2 
(0.56 microns). The color assignments are such that high ratio values are 
red with decreasing values passing through the spectrum ending with the 
lowest values in blue. The saturation of the image is controlled by the 
ratio of TM5 (1.6 microns) over TM7 (2.2 microns). 

TM2 was chosen for its sensitivity to ferric iron oxides; TM7 for its sensi- 
tivity to hydroxyl bands and TM5 for its high variance and broad information 
content. The 5/2 ratio will have high value (red hue) over areas of high 
ferric iron content, vegetation, as well as an assortment of other surface 
materials. The 5/7 ratio will have particularly high values (high saturation 
on the output image) over areas which contain hydroxyl bearing minerals or 
surface materials containing free water (e.g., clays, hydrated salts and 
vegetation). The first eigenband represents a positively weighted sum of the 
seven TM bands and thus provides excellent geomorphologic information allowing 
for precise geographic locations of the image's spectral information. 

We suggest comparison of this image with the 1:250,000 scale Death Valley 
sheet of the Geologic Map series of California. Through comparison with the 
geologic map, some interesting examples of the unique information content of 
the HSV image appear along the northeastern *lank of the Panamint Mountains, 
the eastern Funeral Mountains and the norths n portions of the Renting Spring 
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Mountains. The lower Paleozoic marine section along the northeastern flank 
of the Panamints is clearly distinguished from the older (PC?) section to the 
west. The small outcrops of Tertiary volcanics overlying the Paleozoic 
section are also clearly distinguishable. Note, however, that the Paleozoic 
marine section to the north (Tucki Mountain arec,) is spectrally "confused" 
with the Tertiary volcanics. The Tucki section is distinctly different from 
the Paleozoic sediments to the south of Black Water Wash, however, it is not 
immediately clear why its 5/2 ratio should be so spectrally similar to that 
of the Tertiary volcanics. Along the eastern portions of the Funeral Mountains 
and the Resting Spring Range, there are several examples of st’ .tigraphic 
horizons which are clearly mappable on the HSV imagery and have been grouped 
into the Cambrian marine unit on the 1:250,000 scale geologic map. Although 
such groupings are obviously necessary during geologic mapping, the ability 
to map the individual lithologic beds on the HSV imagery significantly augments 
che information available on the geologic maps. 

SUMMARY 

We in the exploration industry find ourselves in a very challenging situation. 
World consumption of energy and mineral commodities is ever increasing, while 
at the same time, we are at a point where most of the large easy-to-discover, 
cheap-to-produce petroleum accumulation and mineral deposits have been 
located and many of these already exploited and depleted. 

However, on the bright side, our technology is continuing to develop new 
tools with which geologic explorationists can evolve and test new geologic 
concepts. These new concepts allow the geologist to view exploration challenges 
with a new set of glasses, leading to such discoveries as finding oil in 
fractured volcanic rocks in the Great Basin of Nevada and to the thought of 
drilling through igneous and metamorphlc rocks to find underlying oil in the 
Appalachian Mountains. 

The improvements of Thematic Mapper data over multi spectral scanner data 
brings us to the point that we are able to exploit satellite imagery at about 
the same scale that we have used aircraft data in the past. Certainly, 
improved spatial and spectral resolution and wider spectral coverage would be 
welcomed and stereoscopic imagery will be a great boon. However, it appears 
to us that the present TM system offers a near optimum balance between re- 
solution requirements and data handling capabilities. 

In summary, we feel confident that specially enhanced Thematic Mapper imagery 
will make a very significant contribution to the oil and gas and mineral 
exploration industries. The TM's increased spatial resolution will enable 
the production of larger scale imagery, which will greatly increase the 
amount of geomorphic and structural information interpretable. TM's greater 
spectral resolution, combined with the smaller, more homogeneous pixels, 
should enable a far greater confidence in mapping lithologies and detecting 
geobotanical anomalies from space. The results from its applications to 
hydrocarbon and mineral exploration promise to bring the majority of the 
geologic exploration community into that final stage of acceptance and rou- 
tine application of the satellite data. 


Table 1 


COMPARISON OF LAND5AT -1 
SCANNER CHARACTERISTICS WITH 


, -2, -3 AND -4 MULT SPECTRAL 
THOSE OF LANDSAT-4 THEMATIC MAPPER 


MSS 


TM 


Resolution: ^ 80m 


^ 30m 


Wavelength Bands: 




4 a 500 - 600nm 

1 = 

450 - 

520nm 

5 = 600 - 700nm 

2 = 

520 - 

600nm 

6 = 700 - 800nm 

3 = 

630 - 

690nm 

7 = 800 - llOOnm 

4 = 

760 - 

900nm 


5 = 

1550 - 

1 750nm 


6 = 

10400 - 

1 2500nm 


7 = 

2080 - 

2080nm 



IMPLICATIONS OF INFORMATION 
FROM LANDSAT-4 FOR PRIVATE INDUSTRY 


By Dr. John R. Everett and Dr. Jon D. Dykstra 
Earth Satellite Corporation 
7222 47th St., Chevy Chase, MD 20815 


The broader spectral coverage and higher resolution of Landsat- 
4 Thematic Mapper (TM) data open the door for identification 
from space of spectral phenomena associated with mineralization 
and microseepage of hydrocarbon. Digitally enhanced image 
products generated from TM data allow the mapping of many 
major and minor structural features that mark or influence 
emplacement of mineralization and accumulation of hydrocarbons. 
These improvements in capabilities over multispectral scanner 
data will accelerate the acceptance and integration of satellite 
data as a routinely used exploration tool that allows rapid 
examination of large areas in considerable detail. 

(''ere is a fortuitous synchronous development of satellite 
omote sensirg, and the concepts of plate tectonics and verti- 
cal migration of hydrocarbons - two paradigms which are re- 
shaping geological exploration thinking. The steadily improving 
quality of satellite data is providing geologists with a 
valuable tool with which to test, revise and exploit these new 
concepts. 


TM AND GEOLOGIC EXPLORATION 

Since the early 1960‘s, the science of geology has been undergoing a 
major revolution. The new paradigm of plate tectonics and seafloor 
spreading is replacing the older paradigm of a rigid, stable earth. 
Inherent in the acceptance of plate tectonic theory is a growing appre- 
ciation of the role of plate motion in determining the location of 
mineral deposits and hydrocarbon accumulations. It is fortunate that 
developments in spaceborne remote sensing have paralleled these develop- 
ments in geologic thinking. As a consequence, we have remote sensing 
tools that view the earth with appropriate scale and scope to enable us 
to appreciate and map the regional structures that reflect the motions 
of continent-sized segments of the earth’s crust. Wc received our first 
glimpses of the earth from space with photos from the Apollo and Gemini 
flights. The first three Landsat satellites gave us near ubiquitous 
high resolution (80 metre) coverage in four spectral bands. These data 


have had and continue to have enormous impact on all facets of perception 
and management of renewable and non-renewable natural resources and the 
environment. 


Before speculating on the impact of the new data types from Landsat-4, 
it is useful tc take a look at the role data from the first three landsat 
satellites have in geologic exploration and their current level of 
acceptance. In a general sense, Landsat data has made its major contri- 
bution to hydrocarbon exploration in the spatial domain. In mineral 
exploration, Landsat has revealed some spectral Information, but again 
the major contribution is spatial. The synoptic view of over 34,000 
square kilometres of the earth's surface on a single Landsat image 
permits the detection and mapping of major regional structures associated 
with the geologic development of entire geologic provinces. It is also 
possible, through special digital enhancements, to map some of the more 
subtle surface expressions of fracturing, folding and alteration associated 
with buried hydrocarbon accumulations and the emplacement of mineral 
deposits. The data make it possible to interrelate widely separated 
geologic features and detect subtle changes that occur over tens of 
miles and, hence, have gone unnoticed on conventional types of data. 

Perhaps most important of all, the new perspective that the view from 
space provides stimulates us, or forces us, to think of geology in new 
ways and perceive new possibilities. Truly it isn't a panacea but it is 
an extremely powerful tool. It has not by itself "found" an oil field 
or mine, but it has made significant contributions to the exploration 
thinking that led to the discovery of millions of barrels of oil and 
millions of tons of ore. 


At this point, it Is appropriate to consider how industry regards this 
new too 1 . Bob Porter (President, Earth Satellite Corporation) once 
characterized the oil and gas and mineral exploration industries' 
acceptance of Landsat data or any technical innovation as progressing 
through three stages." 

1. An initial "Gee Whiz" stage; 

2. An "Interesting but not yet commercial" stage; and finally 

3. Full acceptance and integration of the innovation as a bona 
fide exploration tool. 

The progression through these stages appears to be a direct function of 
the geologists' and geophysicists' experience and familiarity with the 
tool; in this instance, satellite data. Consequently, not all explora- 
tion groups progress through each stage in unison. 

In the earlier "Gee Whiz" stage, the geologist may be impressed by the 
presence of a familiar structure or landform that is visible on the 
imagery, or by a subtle feature visible on an exotic digitally enhanced 
color image. At this stage, some are tempted to see satellite data as 
the answer to all exploration woes, and the potential to "oversell" its 
capabilities is a real danger. 



This stage is usually quickly replaced by the "Interesting but not yet 
commercial" stage where Imagery produced from the Landsat data is viewed 
mainly as a cheap, low resolution, aerial photo substitute. To so view 
Landsat imagery is to miss both the unique geologic perspective of the 
imagery and the potential contribution available through computer pro- 
cessing and data base integration of the digital data. 

The final stage Is characterized by the integration of the geologic 
information gleaned from the Landsat Imagery with a complete oil and gas 
or minerals exploration model. At this point, the exploration geologist 
fu’ly appreciates the potentials and limitations provided by the satellite 
system and routinely applies the satellite data to the design and exe- 
cution of exploration programs. This whole process takes time. It took 
the gravity meter and the seismograph sbout 40 years to make the transit; 
the digital computer has progressed more rapidly. 

At present, I think a healthy percentage of the exploration community is 
working out of the second stage of acceptance. There is a great promise 
that the Thematic Mapper data will provide the results necessary to 
convince the exploration geologist and, more importantly, the exploration 
financial managers, to move fully into the final stage. However, in 
order for this to occur, there must exist a continuity of data and, 
equally important, an availability of data at a reasonable (that is, 
user justifiable) price. 

The two major advantages of the Thematic Mapper data on Landsat-4 over 
the MSS system are its increased spatial resolution and its greater 
number of narrow, strategically placed, spectral bands (Table 1). These 
new data are so spectacular, they will undoubtedly stimulate a new round 
of the "Gee Whiz" phase. The increased spectral resolution is allowing 
us to map altered zones associated with mineralization based not only on 
iron oxide but on the basis of recognizing rocks and soils rich in 
hydroxyl groups, such as many of the clays formed as a product of hydro- 
thermal alteration. The increased spectral range and resolution promise 
the ability to detect some types of vegetation changes that are associated 
with anomalous mineralization. This will be particularly helpful where 
soil and plants obscure the bedrock. This capability is not definitely 
proven, but theoretically possible and highly anticipated. 

There is a second revolution going on in petroleum exploration geologic 
thinking in addition to the plate tectonic revolution. The old paradigm 
calls for tightly sealed hydrocarbon traps retaining for long periods of 
time (tens or hundreds of millions of years) petroleum that was generated 
and migrated in the distant past. The newly evolving paradigm envisions 
a much more dynamic scenario in which most, if not all, traps leak, and 
the generation and migration of hydrocarbons is a continuing process. 

This implies that there is very little, if any, really old oil or gas, 
rather, only new hydrocarbons generated from old rocks or retained in 
old traps. The hydrocarbon leaked from these imperfect traps moves 
vertically through the overlying rocks to the surface and, in the course 
of its movement, produces a host of chemical changes. The near surface 
environment manifests this leakage in a variety of geochemical, biological. 


geobotanical , or geomorphological anomalies or the simple presence of 
hydrocarbon itself. 


Mineral exploration is not immune to the consequences of this paradigm. 
This peripatetic hydrocarbon has caused the emplacement of a vast amount 
of lead, zinc, uranium, and silver and has potentially played a role In 
localizing some deposits of gold, copper, and barite. This means that 
both petroleum and mineral geologists will need to be cognizant of the 
depositional and thermal histories of areas. 


In the remainder of the talk, I would like to give you an appreciation 
of a few of the ways that the increased spectral and spatial resolution 
of the Thematic Mapper will affect geologic exploration. First, the 
increased spatial resolution. It is clear to us that, with careful 
computer and photographic processing, the quality of the TM digital data 
enables sharp photographic enlargements to a scale of 1:50,000 and, in 
some cases, larger. With clearly interpretable imagery at these scales, 
exploration geologists are able to significantly refine their structural 
interpretations compared to those made from 80 metre resolution MSS 
imagery. 


O 


o 


What was detectable only as a lineamerl on MSS imagery might 
be able to be confidently mapped as a fault. 
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direction of movement along faults will be identifiable in the 
TM imagery, where they are lost in the lower resolution MSS 
Imagery. 


Southern Ontario (Example) 

There are relatively few TM scenes available of areas with strong oil 
and gas interests, however, the July 25 image of Detroit does include 
some of the oil and gas fields of southern Ontario, Canada. 


Thematic Mapper False Color Composite of Southern Ontario 

On the false color image of the Ontario area, we have delineated some of 
the more prominent linears along with the location of the Malden, Col- 
chester and Leamington oil and gas fields. At the Malden and Colchester 
fields, the hydrocarbon accumulations are in fractured, dolomitized, 
Ordovician limestones. The fractures trend WNW. It's a safe bet that 
the lineament marked on the imagery is the surface expression of a major 
through-going structure which is controlling the subsurface fracturing. 

The Leamington field is a little younger and is located in ancient reef 
deposits. Reefs are known to prefer the high edge of structural blocks. 
The intersecting lineaments mapped on the imagery may well define two 
intersecting normal fault zones responsible for the uplift of a block 
edge and the localization of the Silurian reef. 


Cement, Oklahoma (Example) 


The most exciting potential contribution of the TM data is the availability 
of seven carefully placed spectral bands. For oil and gas exploration, 
these spectral bands will be extremely useful for the detection of 
possible surface rock alteration and geobotanical anomalies associated 
with microseepage of hydrocarbons from burled oil and gas accumulations. 

As the leaking hydrocarbons find their way to the surface, they alter 
the chemistry of the rocks through which they pass. At the surface, 
several things can occur: the altered surface chemistry may change the 

spectral and erosional characteristics of the surface rocks and soils 
and/or it may create variations in the density, type or vigor of any 
vegetation growing in the altered soils. 

As a portion of our TM investigations, we are assessing TM's ability to 
identify and delineate areas of surface alteration due to microseepage 
of hydrocarbons at the Cement oil and gas field. The Cement field is 
overlain on the surface by the Rush Springs Sandstone. Rush Springs 
Sandstone is a characteristic red color due to an abundance of ferric 
iron-oxides. The reducing chemical environment associated with the 
hydrocarbon microseepage appears to have altered the insoluble ferric 
iron to soluble ferrous iron, allowing it to be leached out of the 
sandstone. The result is a bleaching of the sandstone overlying the oil 
and gas field. We have only recently received our TM coverage of the 
Cement area. Our spectral analyses are therefore in their earliest 
stages. These next few images are pictures of our interactive image 
processing system and demonstrate our first cut approach at delineating 
the bleached area associated with the Cement field. 

Cement NCC/HSV 

The iron oxide rich areas are delineated on the imagery as areas of 
orange to red color. The quarter frame image shows the location of the 
red coloration following closely the outcrop of the Rush Springs Sand- 
stone. However, the area overlying the Cement oil and gas field is one 
of those areas where the Rush Springs appears to have lost its strong 
iron oxide signature. 

DEATH VALLEY, CALIFORNIA (EXAMPLE) 

The first place to assess potential of TM data to map different rock 
types is in an area of low vegetation density and diverse rock types. 
Within the present range of the TM system there is clearly no more 
vegetation free area than Death Valley, California. 

(Death Valiev NCC/Death Valley FCC) - (Eiqen/HSV) 

The following imagery are of an approximately 1/3 TM scene area of the 
Death Valley, California overpass on 17 November, 1982. The scene 
includes a natural color, false color, eigen and HSV image. 



The Hue, Saturation and Value (HSV) image is one of the more exciting 
images for geologic applications. Through the use of two ratios as hue 
and saturation, and the first elgenband as the value, the resulting HSV 
image possesses the spectral Information of a ratio image and the spatial 
integrity of the first elgenband. 

The hue of the image is controlled by the ratio of TM5 (1.6 microns) 
over TM2 (0.56 microns). The color assignments are such that high ratio 
values are red with decreasing values passing through the spectrum ending 
with the lowest values in blue. The saturation of the image is controlled 
by the ratio of TM5 (1.6 microns) over TM7 (2.2 microns). 

TM2 was chosen for its sensitivity to ferric iron oxides; TM7 for its 
sensitivity to hydroxyl bands and TM5 for its high variance and broad 
information content. The 5/2 ratio will have high value (red hue) over 
areas of high ferric iron content, vegetation, as well as an assortment 
of other surface materials. The 5/7 ratio will have particularly high 
values (high saturation on the output image) over areas which contain 
hydroxy’ bearing minerals or surface materials containing free water 
(e.g., clays, hydrated salts and vegetation). The first eigenband re- 
presents a positively weighted sum of the seven TM bands and thus pro- 
vides excellent geomorphologic information allowing for precise geographic 
locations of the image's spectral Information. 

We suggest comparison of this image with the 1:250,000 scale Death Valley 
sheet of the Geologic Map series of California. Through comparison with 
the geologic map, some interesting examples of the unique information 
content of the HSV image appear along the northeastern flank of the 
Panamint Mountains, the eastern Funeral Mountains and the northern por- 
tions of the Resting Spring Mountains. The lower Paleozoic marine 
section along the northeastern flank of the Panamints is clearly dis- 
tinguished from the older (PC?) section to the west. The small outcrops 
of Tertiary volcanics overlying the Paleozoic section are also clearly 
distinguishable. Note, however, that the Paleozoic marine section to 
the north (Tucki Mountain area) is spectrally "confused" with the 
Tertiary volcanics. The Tucki section is distinctly different from the 
Paleozoic sediments to the south of Black Water Wash, however, it is 
not immediately clear why its 5/2 ratio should be so spectrally similar 
to that of the Tertiary volcanics. Along the eastern portions of the 
Funeral Mountains and the Resting Spring Range, there are several 
examples of stratigraphic horizons which are clearly mappable on the HSV 
imagery and have been grouped into the Cambrian marine unit on the 
1:250,000 scale geologic map. Although such groupings are obviously 
necessary during geologic mapping, the ability to map the individual 
lithologic beds on the HSV imagery significantly augments the information 
available on the geologic maps. 

SUMMARY 

We in the exploration industry find ourselves in a very challenging 
situation. World consumption of energy and mineral commodities is ever 
increasing, while at the same time, we are at a point where most of the 


large easy-to-discover, cheap- toproduce petroleum accumulations and 
mineral deposits have been located and many of these already exploited 
and depleted. 

However, on the bright side, our technology is continuing to develop new 
tools with which geologic explorationists can evolve and test new geo- 
logic concepts. These new concepts allow the geologist to view exploration 
challenges with a new set of glasses, leading to such discoveries as 
finding oil in fractured volcanic rocks in the Great Basin of Nevada and 
to the thought of drilling through igneous and metamorphic rocks to find 
underlying oil in the Appalachian Mountains. 

The improvements of Thematic Mapper data over multi spectral scanner data 
brings us to the point that we arc able to exploit satellite imagery at 
about the same scale that we have used aircraft data in the past. 

Certainly, improved spatial and spectral resolution and wider spectral 
coverage would be welcomed and stereoscopic imagery will be a great 
boon. However, it appears to us that the present TM system offers a 
near optimum balance between resolution requirements and data handling 
capabilities. 

In summary, we feel confident that specially enhanced Thematic Mapper 
imagery will make a very significant contribution to the oil and gas and 
mineral exploration communities. The TM's increased spatial resolution 
will enable the production of larger scale imagery, which will greatly 
increase the amount of geomorphic and structural information interpretable. 
TM's greater spectral resolution, combined with the smaller, more homo- 
geneous pixels, should enable a far greater confidence in mapping litho- 
logies and detecting geobotanical anomalies from space. The results 
from its applications to hydrocarbon and mineral exploration promise to 
bring the majority of the geologic exploration community into that final 
stage of acceptance and routine application of the satellite data. 

Table 1 

COMPARISON OF LANDSAT -1, -2, -3 and -4 MULTISPECTRAL 
SCANNER CHARACTERISTICS WITH THOSE OF LANDSAT-4 THEMATIC MAPPER 


MSS 

Resolution: ^Om 


JM 

^30m 


Wavelength bands: 

4 = 500 - 600nm 1 = 450 - 520nm 

5 = 600 - 700nm 2 = 520 - 600nm 

6 = 700 - 800nm 3 = 630 - 690nm 

7 = 800 -llOOnm 4 = 760 - 900 nm 

5 = 1550 - 1750 nm 

6 = 10,400 - 12,500nm 

7 = 2080 - 2350nm 


SELECTING BAND COMBINATIONS HUH THEM AT I C MAPP.ER DAB- 

by 


Charles Sheffield, 

Ear+h Satellite Corporation, 

7222 47th Street, 

Chevy Chase, Md. 20815. 

Background . Since the human eye employs three primary colors, and the 
Thematic Mapper returns seven bands of data, one obvious problem that arises 
In making color composite Images Is the choice of bands. The choice In 
non-trlvlal, since three bands can be selected from seven In 35 ways. Also, 
any band can be assigned any color. This gives a total of 210 different 
possible color presentations of TM three-band Images. In this note we present 
a way of reducing that „£i0 to a single choice, decided uniquely by the 
statistics of a scene or subscene, and taking full account of any correlations 
that exist between different bands. 

We should remark here that one well-known and widely used approach to 
this problem of choice Is through the use of principal component Images. 
However, such methods offer a new problem as great as the one that they solve. 
For although the first three principal components contain In a statistical 
sense as much Information as can be presented using three colors, the 
resulting scene Is completely data dependent. It Is thus difficult for an 
Interpreter to apply any previous experience of color-surf ace relationships to 
the analysis of a principal components Image. 

Definition of the method . Consider the 7x7 variance-covariance matrix M for 
the scene or subscene. Ignoring for the moment the fact that the thermal band 
Is of Inherently lower resolution than the rest. Any triplet of bands will be 
represented within this 7 x 7 matrix by a 3 x 3 submatrix. 
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Considering now the 3-dimens Jonal subspace spanned by any particular band 
triplet, the associated variance-covariance matrix defines an ellipsoid within 
the subspace. Further, the sum of the squared principal axes of this 
ellipsoid represents the total variance accounted for by these three bands 
(see Figure 1). One could plausibly (but as we shall see, wrongly) argue that 
the best three bands are those with the largest sum of squared principal axes, 
and hence accounting for the largest total variance. This Is, after all, 
exactly the argument applied In employing principal component Images. Since 
the trace of a matrix Is Invariant under rotational transformations, and since 
the sum of squared principal axes Is equal to that trace, the band triplet 
that accounts for the most possible variance can be found from the original 
variance-covariance matrix simply by selecting the three bands with the 
largest diagonal elements. There Is no need to examine all 35 band 
combinations. 

To see what Is wrong with this approach, consider an extreme case where 
there happens to be perfect correlation between a pair of bands. For 
convenience, suppose that those bands are 1 and 2, and suppose that the 
variance of band 1 (and therefore of 2) Is larger than that of any other band. 
The 7x7 matrix M then has the form: 


( a ; 
! a at 



where a>b,c, . . . 



c 



I 


ORIGINAL. PAGE IS 
OF POOR QUALITY 


Sheffield 3. 


The rotation matrix that will diagonalize the upper left 2x2 sibmatrlx 
then has the form: 



and thus after rotation the upper left 2x2 submatrix will have the form: 



As expected, one eigenvalue Is zero; but the other Is the sum of the 'arlances 
from the original bands 1 and 2 . Since a Is assumed to be large, both bands 1 
and 2 will be Included In the triplet that accounts for maximum variance — 
despite the fact that If either one of them Is used, adding the other 
contributes no new Information. 

The problem lies In the use of total variance as the measure for the 
Information content of the band triplets. This Is equivalent to use of the 
sum of squares of ellipsoid principal axes, and there Is no penalty associated 
with a very small principal axis provided that It occurs In association with a 
large axis (see Figures 2 and 3), as was the case for the above example. 

We propose the use of a new measure for the Information content of the 
triplet, and one that avoids the undesirable property demonstrated above. We 
will select the ellipsoid of maximum volume . This discourages selection of 
pairs of bands with high correlation, since In such cases one eigenvalue will 
be close to zero and the corresponding ellipsoid volume will be small. 

Since the ellipsoid volume Is simply 4/3.TTabc, where n, b, and c are the 
principal axes of the ellipsoid, the volume of the ellipsoid associated with a 
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particular band triplet Is a constant multiple of the square root of the 
product of the eigenvalues for the 3x3 variance-covariance matrix of that 
triplet. However, under rotational transformation the product of the 
eigenvalues Is equal to the determinant of the original 3x3 submatrix. Thus 
we can select the band triplet that provides the ellipsoid of maximum volume 
simply by computing and ranking In order the determinants of each 3x3 
principal submatrix of the original matrix M. The bano triplets associated 
with these determinants will then be ranked In order of decreasing overall 
Information content. Given the original matrix M, the total computation to 
achieve this ranking Is trivial. It requires a few hundred multiplications, 
followed by a sort of a list of 35 Items. A BASIC program to perform this Is 
given as an Appendix to this note. 

This procedure gives the best triplet, but the asslgment of colors Is 
still to be made. Now we can make use of the actual variances (the diagonal 
elements of M). Since the eye Is most sensitive to green, next to red, and 
least to blue, we will assign green to the band triplet member of maximum 
variance (I.e. most variation within the Image), red to the triplet member of 
second largest variance, and blue to the triplet member of smallest variance. 
The definition of bands for production of a color Image Is now complete. 

Examples and comments . The procedure has been applied to a number of scenes 
of very different ground cover. Including Washington D.C., Death Valley, and 
Cement, Oklahoma. The results for Washington and for Death Valley are given In 
Tables 1 and 2, together with the associated variance-covariance matrices. The 
following comments apply to all scenes studied to date. 

1) The band combination 1,4,5 (In the order blue, red, green) Is usually, 
but not always, the selected triplet. In cases where It does not rank first, 
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It ranks second or third. 

2) The natural color combination 1,2,3 and the standard false color 
combination 2,3,4, both place far down In the rankings. In the case of 
Washington, the natural color combination Is 29th (lower than anything except 
some thermal band combinations, which are low for another reason to be 
discussed shortly); the 2,3,4 combination was ranked In 16th place. For Death 
Valley, the 1,2,3 natural color combination ranked 32nd, and the 2,3,4 
combination Just above It, at 31st. This Is presumably a consequence of the 
very high correlations between the first four bands. 

3) Triplets that rank high always Include either band 5 or band 6 (note: 
the bands here are ordered by Increasing wavelength, so the thermal band Is 
band 7). This emphasizes the great Importance of these new bands on general 
Information-bearing grounds. 

4) The triplet selected Is not always or even usually the triplet with 
the greatest Individual variances, though large variances are naturally 
preferred somewhat In the selection process. 

Q-ther considerations and comments . 

1) The statistical analysis performed here used P tapes (all that we had 
available) In which the original histograms had already been modified by the 
gains and offsets. It would be preferable to work with data that have had no 
gains or offsets applied, l.e., with A tapes prior to any radiometric 
correction. If band selection of this type becomes common, It would be nice to 
have A tapes generally available from the EROS Data Center. 

2) The thermal band Is of lower resolution than the rest, thus It would 
not be appropriate to give It the same weight In the selection process. How 
should one therefore de-weight It? One argument runs as follows: The maximum 
Information that a scene can contain Is given by the number of pixels, since 
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In the ultimate case there would be no correlation between pixels, and each 
would carry Independent Information about some feature of the surface. In such 
a case, the amount of Information that the thermal band can contribute Is only 
1/1 6th that of the other bands, because there are 16 times fewer pixels In 
that band. Therefore one should deweight the thermal channel by a factor of 
16, Such deweighting was performed In the experiments reported here. However, 
we should also note that this made no difference at all to the preferred band 
triplets, since even without deweighting we found no case where a triplet 
Involving the thermal channel was In the top five. 

3) It Is obvious when one looks at Images created from the triplet 1,4,5 
that for some applications this combination will be much Inferior to others, 
such as natural color and standard false color. This restates the old truth- 
one man's noise Is another man's signal. However, the preferred triplets have 
another advantage: they provide Images of unusual clarity, with far less 
residual striping than Is seen In, for example, the natural color Images. 

4) Although combinations such as 1,4,5 produce Images that are at first 
sight unfamiliar and unusual, the assigned colors are not scene-dependent. 

Thus In contrast to the scene dependent colors of principal component or ratio 
Images, the Interpreter quickly learns to associate colors with particular 
ground condition. We therefore believe that there are definite advantages to 
seeking color composites from the original bands, rather than through band 
ratios or band combinations. 
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Figure Is 


The variance-covariance ellipsoid, 


principal axoa 
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Figure 2. High correlation, 
bands 1 and 2. 


Figure 3 . Low correlation, bands 
1 and 2, but lower 
individual variances. 


2 , 



A, + large, 

*A| . small 



^ +^ 2 . smaller than in Fig. 2. 
• 'A ^larger than in Fig. 2. 


Based on ellipsoid volumes, Fig. 2 case is preferred over Fxg. 1 case,, 
although the former accounts for a greater total variance. 


APPEl'DIX: THE BEST-BAND PROGRAM. 
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20 PRINT "SELECTION OF BEST THREE BANDS BASED ON ELLIPSOID VOLUME" 

30 PRINT "DEATH VALLEY WITH REDUCED THERMAL VARIANCE" 

40 DIM R(36),Q(36) 

50 DIM U( 36) , V( 36) 

60 DIM M(8, 8) 

70 REMARK: M is the varianoe-covarianoe matrix for the scene or subscene. 

80 REMARK: The arrays R,Q,U and V are storage arrays used in the program. 

90 REMARK : Note that the program assumes that band 7 is the thermal data , 
and band 6 is the 2 . 2 micrometre data . 

100 REMARK: The instructions 190 to 230 (except for 220, which sets a count) 
reduce the variance of the thermal channel to allow for the lower spatial 
resolu 1 ' j.j of the thermal channel pixel. 

190 FOR I = 1 TO 6 
200 M(I,7) = M(I,7) / 4 
210 NEXT 
220 C r 1 

230 M(7,7) = M(7,7) / 16 

240 PRINT "RANK DETERMINANT COMBINATION" 

250 FOR I = 1 TO 5 
260 FOR J = I + 1 TO 6 
270 FOR K = J + 1 TO 7 

280 D1 = MCI, I) * (M ( J , J) * M(K,K) - M(J,K) “ 2) 

290 D2 = M(I,J) * (M(J,K) * M(I,K) - M(I,J) * M(K,K) ) 

300 D3 = M(I,K) * (MCI, J) * M(J,K) - M(I,K) * M(J,J)) 

310 DT = D1 + D2 + D3 

315 REMARK: The next instruction makes the determinant an integer; this is 
not necessary, it is done for convenience of output only. 

320 DT = INT (DT) 

330 N = 100 * I + 10 * J + K 

340 R(C) = DT:Q(C) = N 

350 C s C + 1 

360 NEXT 

370 NEXT 

380 NEXT 

385 REMARK: The next piece of code sorts the determinant into descending order. 

390 FOR I s 1 TO 35 
400 N = 0 

410 FOR J = 1 TO 35 
420 IF R(I) | R( J) THEN 440 
430 N = N + 1 
440 NEXT 

450 U(N) = R(I)sV(N) = Q(I) 

460 NEXT 

470 FOR I s 1 TO 35 
480 PRINT I,U(I),V(I) 

490 NEXT 
500 PR# 0 
510 END 
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SELECTION OF BEST THREE BANDS BASED ON ELLIPSOID VOLUME 

THIS IS THE WASHINGTON SCENE WITH REDUCED VARIANCE ON THE THERMAL CHANNEL 


RANK 

DETERMINANT 

COMBINATION 

1 

433858 

145 

2 

205811 

345 

3 

138551 

146 

4 

124784 

245 

5 

101638 

456 

6 

71723 

156 

7 

62960 

346 

8 

49759 

135 

9 

39992 

134 

10 

39609 

246 

11 

36060 

356 

12 

22847 

125 

13 

21953 

256 

14 

16732 

124 

15 

11646 

235 

16 

9709 

234 

17 

7967 

136 

18 

5094 

457 

19 

4752 

157 

20 

3634 

126 

21 

3606 

147 

22 

2294 

467 

23 

2194 

357 

24 

1945 

347 

25 

1616 

236 

26 

1386 

567 

27 

1348 

257 

28 

1130 

247 

29 

727 

123 

30 

688 

167 

31 

276 

367 

32 

215 

137 

33 

175 

267 

34 

84 

127 

35 

43 

237 



SELECTION OF BEST THREE BANDS BASED ON ELLIPSOID VOLUME 
DEATH VALLEY WITH REDUCED THERMAL VARIANCE 


Sheffield 11 



on 


cr\ 

00 

on 

o^ 


on 

in 


on 


OO 

on 

• 

• 

• 

• 

# 

• 

• 

o 

CM 

«=r 

in 

OO 

o\ 

in 

r— 

C\J 

r— 

c» 

•=r 

VO 


on 

VO 

on 

OV 


00 


VO 

i“ 

cn 

o> 

O'* 

on 

on 

• 

■ 

• 

• 

• 

• 

r- 


in 

CM 

00 

VO 

on 

• 

•=r 

o 

VO 

in 

vO 

CM 

CO 


i — ■ 

r* 

*- 

on 

og 





Sheffield 


12 . 


SELECTION OF BEST THREE BANDS BASED ON ELLIPSOID VOLUME 
DEATH VALLEY WITH REDUCED THERMAL VARIANCE 


RANK 

DETERMINANT 

COMBINATION 

1 

1462581 

145 

2 

859695 

156 

3 

684248 

135 

4 

601687 

146 

5 

432952 

345 

6 

346425 

157 

7 

328331 

356 

8 

319827 

245 

9 

275534 

456 

10 

263989 

136 

11 

219239 

256 

12 

204146 

125 

13 

167450 

346 

14 

137060 

357 

15 

127643 

246 

16 

121117 

167 

17 

107494 

457 

18 

103781 

235 

19 

89506 

257 

20 

76827 

126 

21 

75913 

134 

22 

49163 

367 

23 

40621 

467 

24 

39230 

236 

25 

37614 

147 

26 

31621 

267 

27 

21579 

137 

28 

21322 

124 

29 

20256 

567 

30 

9168 

347 

31 

8118 

234 

32 

7895 

123 

33 

7197 

247 

34 

5037 

127 

35 

2407 

237 


Table 2.b 


Ranked results for Death Valley scene. 



